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A striking shape was recently observed for the cellular organelle endoplasmic reticulum consisting 
of stacked sheets connected by helical ramps [T]. This shape is interesting both for its biological 
function, to synthesize proteins using an increased surface area for ribosome factories, and its geo¬ 
metric properties that may be insensitive to details of the microscopic interactions. In the present 
work, we find very similar shapes in our molecular dynamics simulations of the nuclear pasta phases 
of dense nuclear matter that are expected deep in the crust of neutron stars. There are dramatic 
differences between nuclear pasta and terrestrial cell biology. Nuclear pasta is 14 orders of magni¬ 
tude denser than the aqueous environs of the cell nucleus and involves strong interactions between 
protons and neutrons, while cellular scale biology is dominated by the entropy of water and com¬ 
plex assemblies of biomolecules. Nonetheless the very similar geometry suggests both systems may 
have similar coarse-grained dynamics and that the shapes are indeed determined by geometrical 
considerations, independent of microscopic details. Many of our simulations self-assemble into flat 
sheets connected by helical ramps. These ramps may impact the thermal and electrical conductiv¬ 
ities, viscosity, shear modulus, and breaking strain of neutron star crust. The interaction we use, 
with Coulomb frustration, may provide a simple model system that reproduces many biologically 
important shapes. 

PACS numbers: 26.60.-c,02.70.Ns, 87.10.Tf, 87.15.ap, 87.16.D-,87.16.Tb 


Membrane-bound cellular organelles have characteris¬ 
tic shapes, with the endoplasmic reticulum (ER) being 
particularly striking. The ER is an extensive organelle 
displaying three distinct, yet connected, morphologies: 
tubes, sheets and the spherical envelope around the cell 
nucleus. Recent advances in serial sectioning and electron 
microscopy have revealed the stacked ER sheets to be 
connected by helical structures [T]. Just as spiral ramps 
connect the levels of a multilevel parking garage, these 
so-called “Terasaki ramps” allow sheets to connect yet re¬ 
main parallel over scales large relative to the membrane 
thickness. This parking-garage shape is interesting for 
both its biological function, to synthesize proteins using 
an increased surface area for ribosome factories, and its 
mathematical property as the minimizer of a geometrical 
Hamiltonian largely insensitive to details of the micro¬ 
scopic interactions [2] . In these models of membrane me¬ 
chanics, the Terasaki ramps are stable, topological struc¬ 
tures, akin to screw dislocations, that affect the entire 
morphology of an organelle critical to the metabolism of 
eukaryotic cells. 

Independent of considerations from terrestrial biology, 
in this Letter we find very similar shapes arising in our 
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molecular dynamics (MD) simulations of dense nuclear 
phases. These phases are expected in neutron-star crusts 
and in core-collapse supernovae [5]. If matter is com¬ 
pressed to near nuclear densities, competition between 
short-range nuclear attraction and long-range Coulomb 
repulsion can lead to a variety of complex shapes, so- 
called “nuclear pasta” HIS]. Note that Groenewold et 
al. [5] discuss colloidal cluster phases arising from short 
range attraction and coulomb repulsion. For a particu¬ 
lar density and proton fraction, we also find flat sheets 
connected by spiral ramps. In a neutron star, electrons 
scattering from the spiral ramps reduce both the thermal 
and electrical conductivity. This may impact X-ray ob¬ 
servations of crust cooling in transiently accreting stars 
[3] , and may also lead to the decay of neutron-star mag¬ 
netic fields after about a million years [7]. 

There are, to be sure, dramatic differences between 
nuclear pasta and terrestrial cell biology. Nuclear pasta 
has a density near 10^^ g/cm^, fully fourteen orders of 
magnitude more dense than the aqueous environs of the 
cell nucleus. Furthermore, nuclear pasta involves strong 
interactions between neutrons and protons in addition 
to electromagnetic interactions, while cellular-scale bi¬ 
ology is highly-screened, highly-overdamped and domi¬ 
nated by the entropy of water and complex assemblies 
of biomolecules. Nonetheless the strikingly similar ge¬ 
ometry suggests both systems may have similar coarse¬ 
grained dynamics. 

In nuclear physics, the semi empirical mass formula 
predicts the binding energy BE of a nucleus with A nu- 
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cleons (neutrons plus protons) and Z protons [HIH] , 

BE = ayA- - acZ^IA^^^... (1) 

Here a„, a^, Uc are constants describing volume, surface, 
and Coulomb energies. In addition there are other con¬ 
tributions from the symmetry energy and pairing that 
will not be important here. Competition between surface 
and Coulomb energy contributions can lead to complex 
shapes. We perform simulations at a density of n = 0.05 
fm“^. This corresponds to a packing fraction of 5/16 of 
nuclear saturation density, uq = 0.16 nucleons per fm^. 
Here the system forms flat sheets (lasagne) that are con¬ 
siderably thicker than the size of a single nucleon. This 
thickness is determined from a balance of surface and 
Coulomb energies. Note that there are a variety of dif¬ 
ferent shapes with almost the same energy. For example, 
if the density is decreased somewhat, the system forms 
rods (spaghetti) instead of flat sheets and at still lower 
densities the system forms spheres representing isolated 
nuclei. 

Fluctuations about these simple shapes could have 
low excitation energies that may depend on sub-domaint 
terms in Eq. Reinhard et al. m, see also [T3], use a 
leptodermous (thin-skinned) expansion and density func¬ 
tionals in order to calculate a curvature energy term for 
Eq. Mthat goes like QcurveA^^^. The impact of such an 
on nuclear pasta shapes was considered in ref. 
[14j . The curvature energy could be important for nu¬ 
clear fission where a nucleus dramatically changes shape. 
However, Reinhard et al. find they need to calculate the 
energy of very large nuclei with A in the thousands in 
order to extract acurve theoretically. Therefore, this co¬ 
efficient may not be reliably determined from measured 
nuclear binding energies that are known only for a limited 
range of A. 

Instead, we discuss a very different approach that 
has been employed previously in biophysics, but not 
(to our knowledege) in nuclear physics. We consider a 
Helfrich-Canham Hamiltonian Hq nnun] that involves 
a quadratic form in the surface principal curvatures Ci 
and C 2 including both the mean curvature (Ci + C 2 )/‘ 2 , 
and Gaussian curvature C 1 C 2 , 

Ho = ^B J dS{Ci + C2f + B J dSCiC2 . (2) 

Here dS is an integral over the surface area and B 
and B are constants representing effective rigidity mod¬ 
uli. This energy functional was applied in refs. mm to 
biological systems. Our justification for also discussing 
nuclear pasta shapes with the help of Eq. [^is determined 
only after the fact. This equation predicts spiral ramp 
shapes, and their arrangements, very similar to what we 
find in our MD simulations, see below. 

One relatively low energy solution for Eq. [^is lasagne 
with flat surfaces where Ci = C 2 = 0. This gives Hq — 0. 
Another solution involves spiral ramps with Ci = —(72. 
Here the mean curvature is still zero, but the Gaussian 


curvature is negative so that Hq < 0. Thus the Gaussian 
curvature term may stabilize spiral ramp configurations 
in both biological membranes and nuclear pasta. This 
term may also stabilize configurations that have addi¬ 
tional holes such as the “nuclear waffle” shapes found 
in ref. [15]. These shapes consist of flat sheets with a 
two dimensional array of holes. Similarly, the ER dis¬ 
plays a waffle-type pattern which has been referred to 
in the biological literature as “fenestrations” [inj. Al¬ 
ternatively there may be torus or donut-shaped super¬ 
heavy nuclei where the donut shape both reduces 
the large Coulomb energy and is further stabilized by 
the Gaussian-curvature term. 

In this paper we study the self-assembly of these spi¬ 
ral ramp configurations with molecular dynamics simula¬ 
tions of a simple (semi)classical model of nuclear matter. 
Our simulations are for nuclear pasta but they may also 
have implications for phospholipid bilayer membranes. 

Our MD formalism is the same as that used by 
Horowitz et al. in previous works [21 [IHH121 and is briefly 
reviewed here. It is very similar to a model used by others 
[231124j . Our simulation volume is a cubic box with peri¬ 
odic boundary conditions which contains point-like pro¬ 
tons and neutrons with mass M = 939 MeV. Electrons 
are assumed to form a degenerate relativistic Fermi gas 
and are not explicitly included in the simulations. Pro¬ 
tons and neutrons interact via the two-body potentials: 

Vnp{r) = ae~'' + be~^ (3a) 

Vnn{r) = -k (3b) 

Vppir) = (3c) 

r 

The n and p indices indicate whether the potential de¬ 
scribes a neutron-proton, a neutron-neutron, or a proton- 
proton interaction. Meanwhile, r is the separation be¬ 
tween each pair of interacting nucleons and quantities 
a = 110 MeV, b = -50 MeV, c = -2 MeV, and A = 1.25 
fm^ are parameters of the model. Their values were cho¬ 
sen in ref. m to approximately reproduce some bulk 
properties of pure neutron matter and symmetric nuclear 
matter, as well as the binding energies of selected nuclei. 
The screening length A is chosen to be 10 fm. Equation 
describes an intermediate range attraction between p 
and n which binds nuclei and then a short range repul¬ 
sion that causes nuclear matter to saturate at a density 
no = 0.16 fm“^. Finally there is a long range Goulomb 
repulsion between protons. All of the simulations in this 
paper are at a density of n = 0.05 fm“^, a composition of 
40% protons, 60% neutrons, a fixed temperature kT = 1 
MeV, and use a MD time step of 2 fm/c. 

To study self-assembly of these ramps, we start from 
initial conditions where the particle positions are uni¬ 
formly distributed in the simulation volume with a Boltz¬ 
mann velocity distribution. Figure[2shows the configura¬ 
tion of a 40,000 particle simulation at times from 40,000 
to 2,000,000 fm/c. The proton density Up is shown where 
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FIG. 1: (Color online) Self-assembly of a parking garage struc¬ 
ture in a MD simulation with 40,000 nucleons that started 
from uniform random initial positions. Shown in panels (a) 
through (f) are the configuration after simulation times of 
40,000, 400,000, 800,000, 1,200,000, 1,600,000 and 2,000,000 
fm/c respectively. The color map at right is discussed in the 
text. Figure produced with Paraview |25 |. 


each proton is assumed to have a Gaussian shape of 1 
fm^ volume m and the opacity of the colorscale is 0 for 
Up = 0.00 to 0.02 fm“^ and then increases linearly to 
1 at Up = 0.04 fm“^. A light cream color corresponds 
to high density sheets, while lower density surfaces are 
shown in brown. This system undergoes the following 
self-assembly steps: (1) the low density system collapses 
locally to form higher density filaments that meet injunc¬ 
tions, see Fig. (a). This includes the formation of a 
number of topological holes. (2) Next, the filaments start 
to grow to form curved sheets, Fig. [T] (b-c), (3) these 
sheets then start to straighten out over longer length 
scales, Fig. ^ (d). (4) Boundaries between “domains” 
of sheets with different orientation form four left-handed 
and four right handed helical ramps, see Table |Tj The 
sheets straighten out over the full simulation volume, see 
Fig. (e-f), and the ramps move together to form the 
dipole pattern shown in Fig. (a). This pattern has 
four left-handed ramps to the left and four right-handed 
ramps to the right. For this dipole pattern, the ramps 
are seen in Fig. (c) to make about a 45 degree angle 
with the flat sheets. We find this hnal configuration to 
be stable for times of at least 10,000,000 fm/c. 

To study the dependence on boundary conditions, we 
perform simulations with different numbers of particles 
and correspondingly different sized simulation volumes, 
see Table |Ij The smallest simulation, with only 20,000 
particles, forms uniform flat sheets without any spiral 
ramps (not shown). A 50,000 particle simulation forms 
four left handed and four right handed ramps as shown 
in Fig. [^(e). This is also a dipole pattern with the left- 
handed ramps to the left and the right-handed ramps to 
the right. As a result the ramps make a 45 degree angle 
with the sheets, see Fig. [^(f). Finally a 75,000 parti¬ 
cle simulation forms four ramps in a quadrupole pattern 
where one left-handed and one right-handed ramp are to 



FIG. 2: (Color online) Three axis views, X (left), Y(center), 
and Z (right), of the final configuration of three MD simula¬ 
tions. Panels (a)-(c) at top are for the 40,000 nucleon sim¬ 
ulation shown in Fig. 1. Panels (d)-(f) at center are for a 
simulation with 50,000 nucleons and panels (g)-(i) at bottom 
are for a simulation with 75,000 nucleons. The handedness of 
the helical holes in panels (a), (e), and (h) are indicated with 
L for left-handed and R for right-handed. 


the left as shown in Fig. (h). For a quadruple pattern 
the ramps are observed to make a 90 degree angle with 
the sheets, see Fig. (i) and Table [Tj 

The arrangement of the helical ramps shown in Fig. 
[^panels (a), (e) and (h) agrees very well with theoret¬ 
ical predictions in refs. m [2] that are based on Eq. [| 
Guven et al. argue that tension in the sheets leads to 
an effective long range attraction between two ramps of 
opposite chirality that draws them together until a short 
range repulsive bending force stabilizes the pair of ramps 
at a characteristic distance. Some features of our simula¬ 
tions that are in agreement with Guven et al. are: (1) we 
find an equal number of left-handed (L) and right-handed 
(R) ramps, (2) the ramps are all relatively close together 
with a similar characteristic spacing, and (3) ramps in a 
quadrupole pattern make a 90 degree angle with respect 
to the sheets while ramps in a dipole pattern make an 
approximately 45 degree angle with the sheets. These 
common features suggest, after the fact, that Eq. [^may 
also apply to our nuclear pasta model system. 

Most of our simulations form flat sheets connected by 
helical ramps. However we are able to obtain only flat 
sheets, without any ramps, if we add a small one-body 
potential to the system for early times that biases the 
formation of only flat sheets, or start a simulation in a 
smaller box at high densities where the system is nearly 
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TABLE I: Number of ramps Nr in MD simulations with 
different numbers of particles N. 


N 

20,000 

Nr 

0 

Pattern of ramps 

Angle of ramps (deg) 

40,000 

8 

Dipole 

45 

50,000 

8 

Dipole 

45 

75,000 

4 

Quadrupole 

90 


uniform and then very slowly expand the box during the 
simulation until the system reaches the same final density 
of n = 0.05 fm“^. 



FIG. 3: (A) Scanning electron micrograph of a thin slice of 
sheet-like ER from a mouse salivary gland. Serial section¬ 
ing allows three-dimensional reconstruction from large num¬ 
bers of cross-sections (scale bar = 200 nm). (B) 3D recon¬ 
struction of a left-handed Terasaki ramp that appears in the 
black-outlined region in (A). (C) 3D reconstruction of a right- 
handed Terasaki ramp that appears in the white-outlined re¬ 
gion in (A). (Figure adapted from Terasaki et al. [1]) 


We compare our results to biological observations. 
Reconstructing the three-dimensional geometry of ER 
sheets from two-dimensional serial sections reveal that 
the continuity of parallel sheets comes about through 
the helical winding of the “exposed” sheet edge through 
space, i.e. the Terasaki ramp, see Fig. . The core of 
the ramp (cytosolic side of the membrane) has a highly 
negative Gaussian curvature, but potentially small mean 
curvature, given the opposite signs of the two principal 
radii of curvature. No considerations seem to set a pre¬ 
ferred handedness to the ramps, and the observations 
at hand, though statistically small, are consistent with 
right- and left-handed ramps in equal numbers through¬ 
out the organelle. If the sheet edges are treated as ef¬ 
fectively one-dimensional defects (thought to be stabi¬ 
lized by membrane proteins on the cytosolic side), then 
it is sufficient to treat the rest of the membrane by the 
Helfrich-Canham Hamiltonian. A class of solutions min¬ 
imizing that functional are minimal surfaces, which have 
zero mean curvature and locally minimize area. One con¬ 
sequence of the analysis in [2] is that whereas a single 
Terasaki ramp has a logarithmically diverging energy, a 
left-right dipole pair has a finite energy, and a double 
pair of left-right-left-right ramps (a quadrupole) mini¬ 
mizes it further. The existence of ramp dipoles is natural 
in the biological context, though convincing evidence for 
tightly correlated left-right pairs, or for pairs of pairs, is 
currently lacking. As previously discussed, we find both 
left-right pairs and pairs of pairs in our MD simulations 


of nuclear pasta. 

In conclusion, we have performed molecular dynam¬ 
ics simulations using a simple classical model of nuclear 
pasta. Many of our systems spontaneously self-assemble 
to form flat sheets connected by helical ramps. This 
geometry is very similar to that observed in the three- 
dimensional structure of the sheet-like endoplasmic retic¬ 
ulum. Seeing the same helical shapes in the extraordi¬ 
narily different systems of nuclear pasta in neutron stars, 
and in the membranes of eukaryotic cells strongly sug¬ 
gests that these shapes follow from common geometric 
considerations and may be independent of details of the 
microscopic interactions. 

The curvature energy could help stabilize these heli¬ 
cal ramps, suggesting that they may be common in nu¬ 
clear pasta. If true, electron scattering from the ramps 
could reduce the electrical and thermal conductivity of 
neutron-star crust. This may impact crust cooling [5] 
and possibly lead to the decay of magnetic fields [7] . Al¬ 
ternatively the helical ramps, that connect flat sheets, 
may allow paired protons to percolate throughout the 
system and make it superconducting. In addition, the 
ramps may restrict how the sheets can move past each 
other. This could lead to complex viscoelastic behavior 
for nuclear pasta that could dampen r-mode oscillations 
in rapidly rotating neutron stars, see for example |26) . 
Furthermore, the helical ramps likely will increase the 
shear modulus of neutron-star crust and the frequencies 
of crust oscillation modes. These modes may have been 
observed as quasi periodic oscillations during magnetar 
giant flares m- Finally, the ramps likely increase the 
breaking strain, or mechanical strength, of neutron star 
crust. As a result the crust can support larger moun¬ 
tains that, on a rapidly rotating neutron star, efficiently 
radiate gravitational waves |28j . 

One can study soft condensed-matter analogs of nu¬ 
clear pasta in the laboratory even if direct experiments 
on nuclear pasta are not feasible. This could provide 
unique insights. For example, it may be very difficult 
to predict the actual density of ramps in nuclear pasta 
using only first principle simulations. Instead one can ob¬ 
serve the actual density and pattern of Terasaki ramps 
in analog laboratory systems. 

Our simple interaction has Coulomb frustration with 
short-ranged attraction and long-ranged repulsion. This 
seems to provide a simple model system that reproduces 
many biologically important shapes. One reason is the 
fluid nature of bilayer membranes: interactions between 
phospholipid molecules are determined to a large extent 
by the repulsive term. 

Computational advances have made and will make very 
large-scale MD simulations for such systems “easy” and 
these simulations will likely exhibit very rich varieties of 
shapes and phases. Uncovering similarities between dis¬ 
parate physical systems allows connections to be made 
at the deeper level of symmetry, excitations and the ge¬ 
ometry of topological defects. 
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